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INTRODUCTION 

The promise of filamentary composite materials, whose de- 
velopment may be considered as entering its second generation, 
continues to generate intense interest and applications activ- 
ity, Such interest and activity are wellr- founded, since they 
are based on the possibility of using relatively brittle mate- 
rials with high modulus, high strength, but low density in 
composites with good durability and high tolerance to damage 
and which, when they do fail, do so in a non-catastrophic 
manner. Fiber reinforced composite materials of this kind 
offer substantially improved performance and potentially 
lower costs for aerospace hardware. 

Much progress has been achieved since the initial devel- 
opments in the raid 1960 's. Rather limited applications to 
primary aircraft structure have been made, however, mainly in 
a material-substitution mode on military aircraft, except for 
a few experiments currently underway on large passenger air- 
planes in commercial operation and a few military develop- 
ments which have not seen service use. 

To fulfill the promise of composite materials completely 
requires a strong technology base. NASA and AFOSR recognize 
the present state of the art to be such that to fully exploit 
composites in sophisticated aerospace structures, the tech- 
nology base must be improved. This, in turn, calls for ex- 
panding fundamental knowledge and the means by which it can 
be successfully applied in design and manufacture. 



4 


As technology of composite materials and structures 
moves toward fuller adoption into aerospace structures some 
of the problems of an earlier era are being solved, others 
which seemed important are being put into perspective as 
relatively minor, and still others unanticipated or put aside 
are emerging as of high priority. The purpose of the RPI 
program as funded by NASA and AFOSR has been to develop crit- 
ical advanced technology in the areas of physical properties, 
structural concepts and analysis, manufacturing, reliability 
and life prediction. 

Our approach to accomplishing these goals is through an 
interdisciplinary program, unusual in at least two important 
aspects for a university. First, the nature of the research 
is comprehensive - from fiber and matrix constituent proper- 
ties research, through the integration of constituents into 
composite materials and their characterization, the behavior 
of composites as they are used in generic structural compo- 
nents, their non-destructive and proof testing and, where 
the state of the art will be advanced by doing so, extending 
the research effort into simulated service use so that the 
composite structure's long-term integrity under conditions 
pertinent to such use can be assessed. Inherent in the RPI 
program is the motivation which basic research into the 
structural aspects provides for research at the materials 


level, and vice versa. 
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Second, interactions among faculty contributing to pro- 
gram objectives - which is a group wider than that supported 
under the. project - is on a day to day basis, regardless of 
organizational lines. Program management is largely at the 
working level, and administrative, scientific and technical 
decisions are made, for the most part, independent of con- 
siderations normally associcited with academic departments . 
Involvement of this kind includes - depending on the flow of 
the research - faculty, staff and students from chemistry, 
civil engineering, materials engineering and the department 
of mechanical engineering, aeronautical engineering and 
mechanics . 

Both of these characteristics of the NASA/AFOSR program 
of research in composite materials and structures foster the 
kinds of fundamental advances which are triggered by insights 
into aspects beyond the narrow confines of an individual dis- 
cipline. This is a program characteristic often sought in 
many fields at a university, but seldom achieved. 

Overall program emphasis is on basic, long-term research 
in the following categories: (a) constituent materials, (b) 

composite materials, (c) generic structural elements, (d) 
processing science technology and (e) maintaining long-term 
structural integrity. Emphasis has shifted, and can be ex- 
pected to continue to shift from one time period to another, 
among these areas depending on the states of composite mate- 
rials and structures. Progress in the program will be 
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reported in the following pages under these headings. “Those 
computer software developments are also undertaken which 
both support Rensselaer projects in composite materials and 
structures research in the areas listed above and which also 
represent research with the potential of widely useful re- 
sults in their own right. 

In short, the NASA/AFOSR Composites Aircraft Program is 
a multi-faceted program planned and managed so that scien- 
tists and engineers in a number of pertinent disciplines 
will interact to achieve its goals. Research in the basic 
composition, characteristics and processing science of com- 
posite materials and their constituents is balanced against 
the mechanics, conceptual design, fabrication and testing of 
generic structural elements typical of aerospace vehicles so 
as to encourage the discovery of unusual solutions to present 
and future problems. In the following sections, more detailed 
descriptions of the progress achieved in the various component 
parts of this comprehensive program are presented. 

An evaluative site visit was made to the campus on De- 
cember 2 and 3, 1982, headed by Project Monitor Michael Green- 
field. The site visit team, together with members of the 
Composite Structural Materials Industrial Technical Committee 
who attended this meeting, are listed in Table I-l. Discus- 
sion subsequent to briefings and laboratory tours, during 
which progress in the program during the previous six months 
was addressed, provided considerable guidance for the months 
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TABLE I-l 

P ARTICIPANTS IN THE SITE VISIT OF DECEMBER 2 AND 3, 1982 

NASA/AFOSR Evaluators 
Anthony Amos AFOSR/NA 

Michael Greenfield NASA Headquarters 

Robert Johns NASA Lewis 

Norman Johnston NASA Langley 

Mark Shuart NASA Lemg ley 

Industrial Technical Advisory Committee (ITAC) 

Stanley Harvey Program Manager, Composites 

Boeing Commercial Airplane Company 

I„ Grant Hedrick Presidential Assistant for Corporate 

Technology 

Grumman Aerospace Corporation 

Karl Stenberg Branch Chief - Design Composites 

McDonnell Aircraft Company 
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ahead. Salient among the points which emerged are that 
changes in program direction away from those projects in- 
volving generic design concepts to the more fundamental ma- 
terials and structural behavior are likely, as pursued at 
RPI, to have the more significant impact on the applications 
of interest to NASA/AFOSR at this juncture in advanced com- 
posites developments. 

Accordingly, two new program components were planned 
and set up during the ensuing months , one in the Constituent 
Materials area, the other in Composite Materials. Both new 
projects will begin immediately following the current re- 
porting (and budget) period. Professor Bernhard Wunderlich, 
of the Department of Chemistry, will undertake studies to 
determine the precise glass transition temperatures of sam- 
ples from various points within laminates, which will help 
to establish the uniformity of resin matrices at various 
positions within cured parts, as such may be influenced by 
the thickness, or other dimensions of the laminates in ques- 
tion. Professor Ting-Leung Sham will initiate theoretical 
and experimental investigations of the mechanics of delami- 
nation as it occurs at the free edges of loaded laminates. 

Progress in these two new areas will be covered in the 
next semi-annual report. It is expected that projects in 
generic structural elements, processing science technology 
and maintaining long-term structural integrity will either 
v/ind down or reduce the level of their activities as the new 
studies pick up. 
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PART II 

CONSTITUENT MATERIALS 

II-A TRANSVERSE PROPERTIES OF FIBER CONSTITUENTS IN COMPOSITES 
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II -A TRANSVERSE PROPERTIES OF FIBER CONSTITUENTS IN COMPOSITES 

Se;nior Investigator: R. J. Diefendorf 

1 . Introduction 

The transverse properties of reinforced composites are 
of great interest. One reason is that they may impose limi- 
tations in applications. A second reason is that they may 
hold the key to the realization of the full potential of the 
synergistic behavior of composites. In order for the full 
benefits of composites to be attained, a better understanding 
of the constituent materials is required. In the case of 
the transverse direction, composite properties are a function 
of fiber structure and fiber packing. This implies that only 
average transverse properties can be determined. However, 
only very little engineering data on transverse characteris- 
tics of the fiber constituent has been obtained. Current re- 
search efforts are aimed at collecting this data. 

Direct transverse property measurement of the fiber is 
difficult due to the small fiber diameter (typically 8 microns) . 
Hence, unidirectional laminates were used to supply transverse 
property test specimens. By carefully controlling the cure 
process parameters, the volume fractions of the laminates were 
varied while differences in thickness were kept at a minimum. 

The production of these laminates enabled a plot of property 
versus volume fraction to be constructed , from which a fiber 
value in the transverse direction was extrapolated. 
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2. Status 

The dispersion in volume fraction versus number of fields 
for laminates and theoretical arrays has been investigated as 
a first cut at determining homogeneity. As expected, a poor- 
ly fabricated part (in this case molded) has a greater disper- 
sion than even a laminate made by conventional press curing 
(see Figure II-A-1) . Transverse specimen coefficients of 
thermal expansion (reported previously) were analyzed and 
corrected for the quartz push rod used in the dilatometer 
measuring expansion. Transverse specimen modulus measure- 
ments have been performed recently. The concept of homogene- 
ity has helped explain some of the scatter in the transverse 
data results. 

3. Progress During Report Period 

The data for transverse modulus is presented in Table 
II-A-1 and is plotted in Figures II-A-2 through 6. The data 
for the AS-4 and T-300 fibers appears to have similar char- 
acteristics, and this raises the possibility that they can 
be treated as one data group within the experimental scatter. 
The combined data is plotted in Figure II-A-7, and it is 
apparent that the data sets can be treated as one. The rea- 
son for this is revealed by close scrutiny of the constitu- 
ent fiber. Figure II-A-8 illustrates three of the trans- 
versely isotropic fiber types currently available. Figure 
II-A-8a is that of a "wheel-spoke" structure while Figure 
II-A-8b depicts an onion-like structure. A random fiber 




TABLE II-A-1. UNIDIRECTIONAL COMPOSITE MODULUS RESULTS 



T-300 



AS-4 



Kevlar 


Plate 

(1 - v^)/v^ 

l/VfE^ 

Plate 

(1 - Vfj/Vj 

l/VjE^ 

Plate 

(1 - v^)/v^ 


8-T 

. 8349 

1.7277 

11-T 

.9231 
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1-T 

.5174 

‘3.0602 

9-T 

.7668 

1.6284 

12-T 

.6529 

1.1296 

2-T 

.3141 

2.3723 

16-T 

.6181 

1.3039 

21-T 

.6234 

1.0769 

3-T 

.3228 

2.4950 

20-T 

.5432 

1.0746 

25-T 

.5083 

1.0595 

5-T 

.3298 

2.1176 

24-T 

.7241 

1.3428 

36-T 

.6474 

1.3702 

6-T 

.2658 

2.2898 


30- T 

31- T 
34-T 

1.1186 

.6129 

1.1834 

2.2929 39-T 

1.1208 
1.8953 

Type P Pitch 

.5290 

1. 

0954 

HMS 



Plate 
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' f 2 

Plate 
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.5244 

1.8112 
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19-T 

.6556 

1.5353 

14-T 

.6778 

1.3263 


23-T 

.4837 

1.6358 

18-T 

.6556 

1.5107 


33-T 

.7212 

1.699 

22-T 

.5974 

1.4438 


37-T 

2.7037 

6.1316 

26-T 

.5456 

1.4642 





28-T 

1.7933 

3.3019 





29-T 

.7794 

2.0068 





35-T 

1.3810 

2.6388 
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Figure II-A-3. Plot of Derived Modulus Versus Volume Fraction Parameter for AS-4 Fiber Composites 
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Figure II-A-5 . Plot of Derived Modulus Versus Volirne Fraction Parameter for HMS Fiber Composites 
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Figure II-A-7. Superimposed Modulus Data for AS-4 and T-300 Fibers 
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structure is represented by Figure II-A-8c. From first prin- 
ciples, it is intuitive that structures A and B could yield 
similar results for transverse fiber properties. When the 
fibers are placed in a composite and such factors as packing 
and volume percent resin make significant contributions to 
composite behavior, it becomes easy to understand why it may 
well require a radically different fiber structure to pro- 
duce different transverse property results for a given axial 
modulus. The modulus values for HMS and Type P Pitch also 
appear to be similar in nature. A plot (Figure II-A-9) of 
the combined modulus data indicates that , within the results, 
these two data sets can also be treated as one. 

The results of individual sample CTE runs were corre- 
lated to the samples' respective volume fractions (as deter- 
mined in most cases by a Bausch & Lomb FAS/II image analysis 
system) , '.and the results are plotted in Figures II-A-10 
through 14. As in the case of transverse modulus, the data 
can be separated into pairs consisting of the AS-4 and T-300 
(see Figure II-A-15) and the HMS and Type P Pitch (see Fig- 
ure II-A-16) . The cause for concern in the CTE extrapola- 
tions is the inherent scatter in both the CTE and volume 
fraction measurements. In order to achieve a better estimate 
of the fiber property by anchoring the extrapolation at one 
end, the properties of the neat resin were measured. How- 
ever, due to a Poisson's effect that results from the con- 
struction of the matrix in the longitudinal direction at 
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Figure II-A-10. Plot of CTE Versus Volume Fraction for AS-4 Fiber Composites 
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Figure II-A-12. Plot of CTE Versus Volume Fraction for HMS Fiber Composites 
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Figure II-A-15. Superimposed CTE Data for AS-4 and T-300 Fibers 
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Figure II-A-16. Superimposed CTE Data for HMS and Type P Pitch Fibers 
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low volume fractions, the direct resin expansion is not 
directly useful in anchoring the extrapolation. Using S.cha- 

fii * 

pery's ■ equation for longitudinal composite expansion', 
it was possible to determine the restrained expansion case 
and to use this value in the extrapolation. The results for 
the five anisotropic fiber types for both individual and com- 
bined data sets are summarized in Table II-A-2 . The values 
for GTE of T-300 fibers are in good agreement with those 
of Ishikawa and Van Schooneveld (see Table II-A-3) . Little 
data is available for other fibers but they are judged to be 
reasonable based on the T-300 results. 

As previously discussed, it was possible to anchor the 
extrapolation for the GTE values and thereby provide reliable 
estimates. However, it is not possible to perform a similar 
calculation for the modulus of the resin. Instead, to check 
on the reliability of the modulus estimates, statistical 
analysis and a t-value corresponding to plus or minus one 
standard deviation (about 75% confidence level) , the spread 
in the intercept can be calculated. The results are pre- 
sented in Table II-A-4 . The reason for the high values de- 
termined for 1/INT(-) is that as INT(-) approaches zero, its 
inverse (modulus) blows up disproportionately with respect 
to INT(+) . However, inspection of the 1/INT and 1/INT(+) 
terms indicates the spread is not so great as to make the 


Numbers in brackets in this section refer to the references 
which are listed on page 34. 
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TABLE II-A-2 
PAIRED CTE RESULTS 


Fiber 

CTE (E-6) 

I Fiber Pairs 

CTE (E- 

Kevlar 

69.6 

1 

; HMS & P Pitch 

6.8 

HMS 

6.6 

; AS-4 & T-300 

10.0 

T-300 

10.0 



AS-4 

10.3 



P Pitch 

7.8 




SUMMARY OF 

TRANSVERSE FIBER 

PROPERTIES 


Modulus 

CTE 

Fiber 

psi (E+6) 

in/in/C ( 

Kevlar 

.75 

69.6 

HMS / P Pitch 

3.10 

6.8 

AS-4 / T-300 

5.82 

10.0 


TABLE II-A-3 

COMPARISON OF TRANSVERSE T-300 FIBER CTE RESULTS 


CTE 

8.6 X 10“® 
8.8 X 10“® 
10.0 X 10"® 


Temperature 
Range (°C) 


63.9 X 10~® 
53.2 X 10“® 
51.0 X 10~® 


Ishikawa 

Van Schooneveld 

This work 


20-60 

30-98 

35-110 



SCATTER AND INTERCEPT OF EXTRAPOLATED MODULUS RESULTS 


Sample 

Standard 
Error of 
Estimate 

t ^ 70% 
±la 

Standard 

Deviation 

Intercept 

± A 

Int^ 

{+) 

1/Int . 
(-) 

Int . 
(-) 

1/Int. 

(-) 

1/Int. 

Kevlar 

.1582 

.95 

1.33528 

.1503 

1.4856- 

.6731 

1.1850 

.8440 

.7489 

HMS /Pitch 

.3304 

.75 

.32224 

.2478 

.5700 

1.7543 

.0744 

13.4336 

3.1033 

T- 3 00/AS -4 

.1453 

.75 

.17192 

.1090 

.2809 

3.5600 

.0629 

15.9032 

5 . 8166 
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data invalid, and the modulus values obtained are also rea- 
sonable. 


4 . Plans for Upcoming Period 
During the course of this work it became apparent that 
there were a few research areas that might be worth inves- 
tigating. First, the concept of an index of homogeneity re- 
quires more development. Second, Poisson's ratio data, which 
is practically non-existent for the fiber constituent, might 
be a useful piece of information. The effect of nonuniform- 
ities in producing and processing fiber batches on property 
variances might be addressed. Third, the role of the "3rd 
phase" in fiber reinforced composites (i.e., the interface 
between the fiber and matrix) must be more fully developed. 

5. References 

1. Schapery, R. A., "Thermal Expansion Coefficients of 

Composite Materials Based on Energy Principles", Journal 
of Composite Materials, Vol . 2, No. 3, July 1968, p. 380. 
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II I-A FATIGUE IN COMPOSITE MATERIALS 

Senior Investigator; E. Krempl 

1. Introduction 

The deformation and failure behavior of graphite-epoxy 
tubes under biaxial (axial and torsion) loading is being in- 
vestigated. The aim of this research is to provide basic 
understanding and design information on the biaxial response 
of graphite-epoxy composites. 

2 . Status 

Static, axial and torsional strength, elastic moduli 

and fatigue strength under completely reversed load-controlled 

axial loading were measured for graphite-epoxy [±45] thin- 

s 

in * 

walled tubes and reported . A new specimen lay-up [0, 

±45] consisting of four layers oriented in the 45° direc- 
tions sandwiched between an outer and an inner zero-degree 
layer was established. Static tests of this configuration 
were completed. 

3. Progress During Report Period 
Preparation of two technical papers has begun. 

One will cover the static elastic and strength properties of 
the [0, ±45] tubes; the other will be devoted to their fa- 
tigue performance . 


Numbers in brackets in this section refer to the references 
which are listed on page 41. 
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Two batches of six-layer [0, ±45] tubes were prepared 
using the same prepreg and the same curing cycle . The only 
difference in their manufacture was the use of a seal be- 
tween the shrink tube and the aluminum mandrel in the second 
batch. (See Reference [1] for a description of the specimen 
manufacture procedure.) This seal prevented both the bleed- 
ing of epoxy during the curing process and pressure equali- 
zation. 

The results are given in Table III-A-1 and are repro- 
duced as Figure III-A-1. It is seen that each batch has 
different strength properties. The specimens of the first 
batch have generally inferior strength properties when com- 
pared with those of the specimens of the second batch. This 
is surprising, since the only obvious difference in the fab- 
rication procedure is the use of seals for the second batch. 

With the exception of the second quadrant, where low 
strengths are observed for specimens of the second batch (no 
specimens of the first batch were tested in this direction) , 
the failure points seem to form an ellipse. 

Compared to the properties of the [±45] 4-layer tubes 

s 

reported in Figure III-B-1 of Reference [2] , the addition of 
the layers of 0° fibers has roughly doubled the axial tensile 
failure stress and has caused a difference between the ten- 
sile and compressive strengths. The marked asymmetry in the 
torsional strengths shown in Figure III-B-1 of Reference [2] 
for the [±45]^ tubes has disappeared. However, the torsional 
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TABLE III-A-1 


BIAXIAL 

STATIC 

FAILURE 

STRENGTH 

OF [0, ±45] GR/E TUBES 
s 

Tube 

No. 

Stacking 

Sequence 

a 

MPA 

T 

MPA 

Failure Location Eatch 

Upper Central Lower No. 

141 

[0/ 

±45] 

s 

-262 

0 

X 1 

142 

II 

II 

0 

-132 

X " 

143 

II 

II 

0 

+159 

X " 

144 

II 

II 

+290 

0 

X " 

145 

II 

II 

-161 

-93 

X " 

146 

II 

II 

+218 

+127 

X " 

1 

II 

II 

-116 

+65 

X 2 

2 

II 

II 

+ 147 

-15 9 

X " 

3 

II 

II 

0 

-181 

X " 

4 

II 

II 

0 

+130* 

XX " 

5 

II 

II 

+ 341 

0 

X " 

6 

II 

II 

0 

+191 

X " 

7 

II 

II 

-271 

0 

X " 

8 

II 

It 

+ 252 

+144 

X " 

9 

II 

II 

-145 

+ 84 

XXX " 

10 

II 

II 

-221 

-126 

X " 


Preloaded Specimen 
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fciilure stress of the [0, ±45] ^ tubes is smaller than the 

failure torsional stress of the [±45] tubes. (Failure 

s 

stress properties are based on the total, i.e., fibers plus 
resin, original cross-sectional area of the tubes.) 

4. Plans for Upcoming Period 
The static and fatigue properties of [0, ±45]^ tubes 
will be explored further. The two papers mentioned earlier, 
stemming from the Doctor of Engineering thesis of T.-M. Niu, 
will be prepared for submission to an appropriate journal. 

5 . References 

1. Krempl, E. and T.-M. Niu, "Graphitej-Epoxy [±45]s Tubes. 
Their Static Axial and Shear Properties and Their Fatigue 
Behavior Under Completely Reversed Load Controlled Load- 
ing*', Journal of Composite Materials, Vol. 16, 1982, 
pp. 172-187. 

2.. 42nd Semiannual Progress Report, Composite Structural 
Materials, RPI, Sept. 30, 1981 - April 30, 1982, July 
1982. 


6. Current Publications or Presentations by 
Professor Krempl on this Subject 

"Biaxial In-Phase and Out-of-Phase Behcivior of Graphite- 
Epoxy Tubes" 

Presented at the ASTM Conference on Biaxial/ 
Multiaxial Fatigue, San Francisco, California, 
December 15-17, 1982. 
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III-B MATRIX DOMINATED PROPERTIES OF HIGH PERFORMANCE COM- 
POSITES 

Senior Investigator: S. S. Sternstein 

1 . Introduction 

This project is concerned with those properties of high 
performance composites which are strongly dependent on the 
polymeric matrix. 

Moisture is known to adversely affect the properties of 
both neat epoxy resin and epoxy matrix composites. A major 
contribution to the moisture degradation of mechanical prop- 
erties in neat epoxy is due to the internal stresses and 
strains caused by inhomogeneous swelling. The swelling is 
inhomogeneous since neat resin contains structural fluctua- 
tions, or nodular regions of nonuniform crosslink density, 
which are produced during curing. Several factors which 
contribute to the nodular structure include improper mixing 
of the reagents, thermodynamically based partial segrega- 
tion, excessive intramolecular reaction and incipient forma- 
tion of nodules before macrogelation. 

In addition, swelling is inhomogeneous in composite 
laminates because of the additional nonuniformity of the 
system due to the fibers or any other filler having swelling 
properties different from those of the matrix. 

2 . Status 

To date, specific investigations have included visco- 
elastic characterization of the glass transition region in 
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both laminates and neat resins, delamination studies, mois- 
ture interactions and inhomogeneous swelling phenomena. 

3. Progress During Report Period 

During the present report period, our objective has 
been to further quantify both experimentally and theoreti- 
cally the interactions of water with epoxy based systems and 
other amorphous, glassy polymers. 

a. Experimental Studies 

A set of composite samples, exposed for various periods 
of time to a 100 °C, 100% relative humidity environment (as 
opposed to the boiling water environment described in the 
previous semiannual report) , were tested, and the in-phase 
(storage) stiffness (M') results are shown in Figure III-B-1 
versus weight gain (moisture) . The values of M* are norma- 
lized by the value of M' for the dry sample and also norma- 
lized by the cube of the thickness for the wet sample (as 
tested). It is clear that, at about 2.5% moisture pickup, 
the close of the curve markedly increases , suggesting a 
change in the structure-moisture interaction mechanism. The 
loss factor, M"/M' , for composite samples is given in Figure 
III-B-2. The ordinate represents the increase (or increment) 
in loss factor which is obtained in a wet sample relative to 
a dry sample. If moisture content is increased only (no 
drying) , all samples indicate a linear increase of loss 
factor with moisture content. If the sample is dried, then 
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Figure III-B-1 

Normalized In-Phase Stiffness as a Function of Moisture Weight Gain for Samples Exposed to Water Vapor 
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a permanent residual increase in loss factor (tan 6) is ob- 
served. The same behavior was reported previously, when 
composite • samples were exposed to a boiling water environ- 
ment. This indicates that the moisture effects on dynamic 
mechanical properties of composite materials are the same 
regardless of different states of moisture environment. 

The hysteresis effect in loss factor (shown in Figure 
III-B-2 and the previous report) indicates that composite 
samples are permanently damaged. The micrograph of a com- 
posite sample with 5.1% weight gain is shown in Figure III- 
B-3. The damage is easily visible at 50 magnafication. 

Note that cracks occur mostly at interfaces between laminae, 
with a small amount of interply crackings. The micrograph 
suggests that the large weight gain and thickness increase 
are due to moisture trapped in the separated laminae . It 
is clear that the reduction in stiffness is due to the delam- 
ination induced by moisture uptake . 

During the previous period, we reported that the change 
in M' of neat resin occurs at about 7% moisture uptake. 

Since the composite samples are about 35% resin, this sug- 
gested that the slope change in the composite is strongly 
related to the failure of the matrix and is therefore matrix 
controlled. 

In order to further investigate the damage mechanism 
induced by moisture uptake, the polarized light microscope 
was used to provide microscopic evidence of the failure 




Figure III-B-3. Micrograph of Composite with Moisture Uptake of 5.1% 
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patterns. The micrograph at 100 magnification of a dry, thick 
epoxy resin (~ 2 ram) sample when viewed with polarized light, 
is shown in Figure III-B-4. The light intensity pattern sug- 
gests that the epoxy resin is not a homogeneous material. As 
a matter of fact, the pattern is an overlapping image of mul- 
tiple layers of epoxy resin. The pattern is diminished when 
moisture uptake is about 6% as shown in Figure III-B-5. 

Further moisture uptake (about 8%) actually causes the sample 
to produce microcracks as shown in Figure III-B-6. 

A corresponding study in thin films (~ 0.1mm) reveals 
that inhomogeneous swelling occurred in the epoxy resin upon 
exposure to moisture. The 100 magnaf ication of a dry, thin 
film shows a distinct light intensity pattern as shown in 
Figure III-B-7. It is believed that the bright areas are 
the precipitates of curing agent (dycyanodiamide) of the 
cured epoxy. The birefrigence patterns are observed (gray 
areas) due to the thermal mismatch between the precipitate 
and the surrounding matrix (network of epoxy resin) . The 
light intensity pattern produced by swelling in the vicinity 
of precipitate, when the sample is exposed to 20-hour boil- 
ing water environment, is shown in Figure III-B-8. The com- 
plete extinction along the polarizer axis (vertical) and 
analyzer axis (horizontal) , approximate symmetry of the four 
quadrants, maximum interfacial retardation and decay of re- 
tardation with distance from the interface, is typical of 
birefringence patterns which retained their interfacial in- 
tegrity upon swelling. 
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Figure III-B-5. Birefringence Pattern of Epoxy with Moisture Uptake of 6.2% 



Figure iII-B-6. Birefringence Pattern of Epoxy with Moisture Uptake of 7.9% 
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— — 100 i* lOOX 

Figure III-B-7. Birefringence Pattern of Dry Epoxy Film 




100/^ 100X 


Figure III-B-8. Birefringence Pattern of Epoxy Film 

Exposed to Boiling Water for 20 Hours 



The changes in birefringence pattern, when the epoxy 
samples are exposed to a boiling water environment for vari- 
ous periods of time (i.e., 120 hours, 370 hours and 440 
hours), are shown in Figures III-B-9, 10 and 11, respective- 
ly. In Figures III-B-9 and 10, more birefringence patterns 
are observed, and they suggest that the epoxy resin is so 
inhomogeneous that the diffusion and swelling process are 
strongly dependent on the local structure of the material . 
Previously, other investigators have suggested that inhomo- 
geneity of crosslink density in cured epoxy exists, based on 
electron microscopy of fracture surfaces. It has been known 
that, if the crosslink density was insufficient to produce a 
network structure capable of withstanding the high radial 
extension ratio at the interface, a cohesive failure in the 
interfacial region developed. This is shown, in Figure III- 
B-11, by the disappearance and/or relaxation of light inten- 
sity pattern, when epoxy is exposed to boiling water for 440 
hours. The difference in birefringence patterns between 
Figures III-B-8 and 11 indicates that the cohesive failure 
occurred in the inter facial region. 

The dynamic mechanical data and microscopic evidence 
suggests that the effect of moisture and temperature on the 
mechanical performance of composites and neat resin is detri- 
mental, that the matrix stiffness is a significant factor in 
the out-of-plane stiffness of composites and that inhomogeneous 
swelling in the epoxy matrix is the major mechanism responsible 
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Figure III-B-9, Birefringence Pattern of Epoxy Film 

Exposed to Boiling Water for 120 Hours 
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Figure III-B-10. Birefringence Pattern of Epoxy Film 

Exposed to Boiling VJater for 370 Hours 
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Figure III-B-11. Birefringence Pattern of Epoxy Film 

Exposed to Boiling Water for 440 Hours 
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for the degradation of mechanical properties of composijbes 
and neat resins upon exposure to moisture. 

b. Theoretical Studies 

Using the thermodynamic concepts of swelling and mix- 
ing, a constitutive equation relating stress, strain and 
volume fraction of liquid in the swollen polymer has been 
developed. This nonlinear constitutive equation is appli- 
cable in the reversible regions of swelling to amorphous 
network polymers in swelling equilibrium with a solvent at 
constant temperature and constant external pressure. Equa- 
tions for inhomogeneous swelling in polymeric materials with 
a spherical inclusion havej then been developed for the case 
of zero external load by minimizing the total Gibbs' free 
energy. The spherical inclusion, here, may represent a filler 
particle dispersed in a structural fluctuation field. It is 
assumed that the inclusions are isotropic and stress-free 
in the unswollen state and are located far from each other 
so that they do not interact with the inhomogeneous swelling 
field of others. The numerical solution of the equations 
for Narmco 5208 epoxy resin and 100% relative humidity at 
25 °C surrounding environment is given in Figures III-B-1 
through 6 for a few cases of inhomogeneities. 

Figures III-B-12 and 13 show the stresses and volume 
fraction of water, respectively, as functions of reduced 
rcidius (starting from the inclusion interface ) for the case 
of a rigid (nonswelling) spherical inclusion imbedded in and 
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Figure III-B-12. The Radial and Tangential Stresses in a Swollen Epoxy- 
Matrix of Infinite Extent as Functions of Reduced Dis- 
tance From a Rigid Spherical Inclusion. Inset shows 
the schematic of stress field. 



REDUCED RADIUS, r 


FigureIII-B-13 . The Volume Fraction of Water as Functions of Reduced 
Distance for the Same Case as Figure III-B-12. 
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bonded to a homogeneous, infinite epoxy matrix. The reduced 
radius is the radius normalized with respect to the inclu- 
sion radius, R. Stress and volume fraction of water are 
maximum at the inclusion interface, thus indicating that any 
damage would most likely start there. 

Because of limited knowledge about the mechanism of 
formation of inhomogeneities, in neat or in situ resin and a 
lack of experimental work, a precise characterization of the 
structural fluctuation field is difficult. At present a 
single parameter, p, defined as a material property of which 
the bulk and shear moduli of the neat resin are a linear 
function, is assumed to characterize the structural fluctua- 
tion field. Also, it is assumed that the fluctuation field 
is spherically symmetrical and thus can be expressed as a 
continuous or discontinuous function (i.e. , Gaussian, expo- 
nential, step, etc.) of the reduced radius. The reduced 
radius in this case is radius normalized with respect to some 
characteristic radius like the radius where fluctuation 
(|p - l|) first reaches a fixed fraction [here exp (- 2 )] of 
maximum fluctuation. Figure III-B-14 shows the stresses as 
functions of reduced radius (starting from the center of the 
fluctuation field) for a positive Gaussian fluctuation field 
(shown in the inset) . The corresponding distribution for 
volume fraction of water is shown in Figure III-B-15. As 
expected, it is seen that the stress fields produced in neat 
resin by a large structural fluctuation field (large p) of 
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REDUCED RADIUS, r 

Figure III-B-14. The Radial and Tangential Stresses Associated with a 
Positive Gaussian Fluctuation Field in Neat Resin, as 
Functions of Reduced Radius. Inset shows the nature of 
fluctuation field. 



REDUCED RADIUS, r 

Figure III-B-15 . The Voliame Fraction of Water as Functions of Reduced 
Distance for the Same Case as III-B-14. 
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step function profile is the same as the stress field pro- 
duced by a rigid filler particle (Figure III-B-12) . Figure 
III-B-16 shows the stress field for a negative Gaussian fluc- 
tuation field, and Figure III-B-17 shows the stress field 
for a positive fluctuation field dispersed in a negative 
fluctuation field. As is obvious from all these figures, a 
small structural fluctuation field can produce a large stress 
concentration due to swelling, and these stresses alone or in 
conjunction with a stress field due to an external load (The 
two are not simply additive.) can initiate damage at the in- 
terface . 

4. Plans for Upcoming Period 
The experimental moisture interaction studies will be 
continueid using postcured samples. Improvements in the the- 
oretical constitutive equation and a more realistic charac- 
terization of structural fluctuations will be sought. 


5. Current Publications or Presentations by 
Professor Sternstein on this Subject 


"Characterization of the Matrix Glass Transition in Carbon- 
Epoxy Lciminates using the CSD Test Geometry" , with P. Yang. 

Published in The Role of the Polymeric Matrix in 
the Processing and Structural Properties of Com- 
posite Materials, J. C. Saferis and L. Nicolais, 
eds. , 1983 . 

"Dynamic-Mechanical Response of Graphite/Epoxy Composite 
Laminates and Neat Resin", with P. Yang and L. Carlsson. 

Published in Polymer Composites, Vol. 4, No. 2, 

April 1983. 

"Characterization of High Performance Composites" 

Presented Invited Plenary Lecture to the High 
Polymer Physics Division of the American Chem- 
ical Society, Los Angeles, CA, March 22, 1983. 
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Figure III-B-16, The Radial and Tangential Stresses Associated with a 
Negative Gaussian Fluctuation Field in Neat Resin, 
as Functions of Reduced Radius. Inset shows the 
nature of fluctuation. 
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REDUCED RADIUS, r 


Figure III-B-17. The Radial and Tangential Stresses Associated with a 
Fluctuation Field of Nature Shown in the Inset 
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III-C NUMERICAL INVESTIGATION OF MOISTURE EFFECTS 

Senior Investigator: M. S. Shephard 

1 . Introduction 

This analysis builds on the theoretical developments of 
Professor S. Sternstein (as reported in the semiannual pro- 
gress reports on this program, dated December 1979 and June 
1980) , in which the fully-coupled, non-linear thermomechani- 
cal equations of inhomogeneous swelling of composites in the 
presence of moisture and temperature are analyzed. In the 
research reported here, the; one-dimensional case of a single 
fiber in an infinite matrix has been extended to two-dimen- 
sional, multiple fiber cases with far field stresses. 

2 . Status 

As discussed in earlie;r progress reports, moisture ef- 
fects were introduced into the problem through a nonlinear 
constitutive relation, the principle of virtual displacements 
was employed to develop a nonlinear matrix equation using 
displacement-based finite elements, and these nonlinear equa- 
tions were solved using quasi-Newton finite element analysis. 
Subsequently, the constitutive equation was built into the 
finite element analysis. Test cases with multiple fiber 
packings have been analyzed, in the presence of various hu- 
midity le;vels and applied far field stresses. 
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3. Progress During Report Period 
During the current report period, additional cases were 
run using -the subject analyses, which is the work of graduate 
student Frida Lumban-Tobing. The problems run included both 
one and two fiber problems considering additional loadings 
and hiimidities. The results confirm the trends previously 
reported and are discussed in detail in Dr. Lumban-Tobing ' s 

rn* 

thesis ^ ^ . 


4. Plans for Upcoming Period 
These studies are now considered to be complete, and 
the only remaining tasks are to document the computer program 
for future use and modification. 


5. References 

1. Lumban-Tobing, P. E. R. , "Finite Element Analysis of 

Moisture and External Load Effects in Composites", Ph.D. 
Thesis, Rensselaer Polytechnic Institute, Troy, NY, May 
1983. 


6. Current Publications or Presentations by 
Professor Shephard on this Subject 

"Finite Element Analysis of Moisture Effects in Graphite- 
Epoxy Composites", with F. E. R. Lumban-Tobing and S. S. 
Sternstein . 

Published in Computers and Structures, Vol. 16, No. 1-4, 
1983, pp. 457-469. 

Presented at the Symposium on Advances and Trends 
in Structural and Solid Mechanics, Washington, D.C., 
October 4, 1982 (Given by F.E.R. Lumban-Tobing). 


Numbers in brackets in this section refer to the references 
which are listed on page 54. 
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III-D NUMERICAL INVESTIGATION OF THE MICROMECHANICS OF COM- 
POSITE FRACTURE 

Senior Investigator: M. S. Shephard 

1. Introduction 

To understand the mechanisms of failure in composites 
it is necessary to develop insight into their micromechanical 
behavior, including interactions between matrix and fibers 
as the load is increased from zero to that corresponding to 
failure. Investigations of these phenomena, either experi- 
mental or numerical , are difficult. The purposes of this 
project, being carried out by graduate student Nabil Yehia, 
are to develop the nonlinear finite element analysis capa- 
bility required for and to perform numerical investigations 
of the micromechanical failure of composites. 

2. Status 

The overall structure of a general purpose finite ele- 
ment code for fracture problems has been designed (under this 
project) and is currently being implemented and tested. The 
major requirements for the program are that it (a) accommodate 
various types of fracture problems and (b) house a variety of 
algorithmic approaches to fracture and debonding analysis. 

The linear portion of the program has been coded and tested. 
The fracture criteria are currently being implemented and 
tested. The automatic crack-tracking algorithms have been de- 
veloped and are being tested in a stand-alone mode by graduate 
student Gary Burd. 
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3. Progress During Report Period 
During the last reporting period, effort concentrated 
on two areas of development; implementation and testing of 
numerical algorithms required for the fracture problem and 
development of a meshing algorithm to automatically generate 
finite models accounting for crack growth. Efforts on the 
fracture problem have focused on procedures to predict crack 
growth increments and to insure the level of accuracy ob- 
tained in the calculation of crack-tip stress intensity 
values. The meshing algorithm being developed to automatic- 
ally track crack growth employs the modif ied-quadtree mesh 
generator discussed in the semi-annual progress report for 
this program, dated July 1982, in an algorithm that automa- 
tically meshes a geometry after it has changed due to crack 
growth . 

a. Stress Intensity Factors Calculation 

The stress intensity factors, and , are evaluated 

from the local displacements of the two crack-tip/crack- 

surface singular elements , using the displacement correla- 
[11 * 

tion technxque . The proper elastic, crack-tip singular 
element can be obtained from the standard quadratic isopara- 
metric element by appropriate positioning of the node 
[2 3] 

points ' . Tests have been run on both the singular qua- 

dratic isoparametric triangle (SLST) and collapsed quadratic 
isoparametric quadrilateral (SQ8) . 

* 

Niambers in brackets in this section refer to the references 
which are listed on page 82. 
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To test this module, three sample problems for fracture 

mode 1, as suggested in Reference [4], are analyzed. The 

geometry and dimensions used in the analysis are shown in 

Figure III-D-1. The portions of the geometry analyzed and 

conditions applied for the double edge-notched plate are 

shown in Figure III-D-2. Figure III-D-3a shows a typical 

element mesh used to analyze a problem while Figure III-D-3b 

shows a typical mesh at a crack tip. Although it has been 

[1 2 ] 

repeatedly reported in numerous meetings and papers ' 
that the element side-length of the crack tip affects the 
accuracy with which the crack-tip stress intensity factors 
can be calculated, we have found in this research that this 
effect is almost within the level of engineering accuracy 
for the range of element side lengths under half of the mesh 
length. Figure III-D-4 shows the percent error in evaluating 
for different element side lengths, for the three sample 
problems . It has been also found that the strain energy 
density field is not a strong function of the element side 
length for a reasonable range of edge lengths , as shown in 
Figures III-D-5 and 6. 

b. Crack Propagation Module 

The purpose of this module is to determine the direction 
and magnitude of incremental crack growth, after it has been 
determined that the crack will propagate. Two criteria have 
been chosen as candidates for this module; namely, the minimum 
strain energy density criterion and the maximiim dilatational 
strain energy criterion. 











Figure III-D-3. 


Typical Finite Element Mesh 
b. Closeup of Mesh at Crack Tip 


H' 




-3 


= 1.007 ksi(in)^ 

!L = Length of singular 
element edge 


Figure III-D-4. Error in Computing as a Function of Element Side Length 
b. Central notched plate 






Figure III-D-5. Strain Energy Density Field for Element Side Length Equal To 20% of Cracked Length 





strain Energy Density (dw/dv) . (10)^* ksi 
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The minimum strain eneirgy density has been used^^^ 

to predict the crack propagation increment. This approach 

assumes that the critical strain energy of the material is 

a matericil property and is known for the material under 

study. The latter assumption is not easy to live with in 

most cases, since it requires experimental determination of 

the strain energy density at fracture for each material. A 

[71 

procedure has been suggested for determining this value. 

Once the material critical strain energy is known, the crack 

propagation angle and the crack-tip stress intensity factors 

are determined; the propagation increment can be determined 

as shown schematically in Figure III-D-7. 

The maximum dilatational strain energy criterion has 
[8 91 

been suggested ^ to determine the angle of crack propaga- 

tion. This criterion defines a material critical parameter, 
"^vcr' particular value of the proper point as the 

value of this quantity which indicates that the crack will 
propagate;. The Von Mises elastic-plastic boundary is used 
as the curve on which the maximum value of T is evaluated. 

V 

To determine the crack propagation angle , one can decompose 
the: total strain energy at a point a distance, r, from the 
crack tip as follows: 

T = T + T^ 

V d 

T^ = the dilational strain energy component 
T^ = distorsional strain energy component 


where 



strain Energy Density (dw/dv) . (10)** ksi 





79 


and each can be defined as follows: 
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T 


V _ 

ifb, .K^ 
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r 

r [ 11 I 
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where : 

Ki , and are the crack-tip stress 

intensity factors and. 

b. . and c. . are functions of the local 
i: i: 

angle 0 and the material properties. 

Although this criterion has been used just to determine 
the crack propagation angle it can be extended to de- 
termine the magnitude of the crack propagation increment as 
well. To determine the crack propagation increment, one can 
use either the critical dilatational strain energy approach 
or the critical distance to the elastic plastic boundary and 
check against the appropriate value of the predicted direction 
of propagation on the elastic-plastic boundary around the 
crack tip. Unlike the minimum strain energy density criter- 
ion ^ ' ,, both critical values mentioned above can be easily 

determine;d for a given material provided that the material 
fracture toughness is given. 

The T criterion has been implemented and preliminarily 
tested on the sample problems mentioned earlier in this re- 
port. Reasonable results were obtained. Further testing is 


required , however . 
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c. Error in Predicting the Crack Increment 

In Figures III-D-5 and 6, the error in is defined as; 


K - K. 


(K ) 

I ER 


IC 


K 


IC 


where : 


is the calculated value and 
is the critical value. 

The error in predicting the crack propagation increment, 
using either of the criteria discussed above, can be obtained 
as follows : 
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This implies that the error in predicting the magnitude 
of the crack propagation increment is about twice as much 
as evaluating the crack-tip stress intensity factor. In 
other words, the shorter the singular element side length, 
the better will be the prediction of the crack increment, 
based on the results obtained > with a range of L/a between 
0.0 and 0.3 required to yield results within ten percent 
(reasonable engineering accuracy) . 


d. Automatic Mesh Generation for Crack Tracking 

The ability to automatically track the crack growth re- 
quires an automatic mesh operator and an adaptive 

geometry data base. Alterations are being made to the modi- 
f ied-quadtree automatic mesh generator so that it can be used 
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for this purpose. Thus far, work has been done in three 
major areas. These areas are (a) converting the input of 
the mesh generator so as to work from an adaptive geometry 
data base, (b) allowing a quadrant to have several boundar- 
ies within it and (c) adding pointers to the quadtree. 

Because several changes of geometry will have to be 
made during the analysis of a given problem, an adaptive 
data structure is required. Linked lists lend themselves 
quite nicely to this job. A region can be represented as an 
area enclosed by a linked list of curves; and a curve can be 
represented by a linked list of x-y coordinate pairs. As 
the geome'.try changes , only a few pointers need to be changed 
rather than reorganizing the storage of data in memory. The 
curve discretization algorithm, used by the present version 
of the modif ied-quadtree program, has been recoded to work 
from this linked-list data structure. 

The original version of the modif ied-quadtree program 
allows only one boundary to be represented within a quadrant. 
At least two boundaries will always have to be represented 
in the analysis of crack .growth. In order to make this task 
more manageable, the possible number of cutting points in a 
quadrant has been limited to four. If the diagonals used to 
convert quadrilaterals into triangles are chosen properly, 
the limitation on the number of cutting points does not ad- 
versely affect the quality of triangular meshes. An algorithm 
for choosing diagonals has been successfully implemented. 
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Because several meshes will be generated during a crack 
growth problem, it is important that the mesh generation pro 
cess be as inexpensive as possible. A significant improve- 
ment in the cost of generating a mesh can be made by keeping 
pointers from a parent to its sons , rather than by deriving 
this information by a tree traversal every time it is needed 
Accordingly, portions of the mesh generator are being re- 
written so that it will use pointers. 


4. Plans for Upcoming Period 

During the next reporting period, effort will concen- 
trate on the following: 

a) interfacing the automatic mesh generator for 
remeshing as cracks propagate with the analy- 
sis routines, 

b) continuing the testing of crack propagation 
modules and 

c) carrying out specific micromechanical analy- 
sis on composites that can be properly repre- 
sented with those features available in the 
program. 
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III-E ADVANCED ANALYSIS METHODS 

Senior Investigator: E. J. Brunelle 

Introduction 

The effort to find solutions for individual lamina in 
terms of a minimum number of generalized variables has 
gradually widened to finding solutions for composite plates 
with commonly used stacking sequences. In particular, anti- 
symmetric cross-ply and anti-symmetric angle-ply plates 
(for particular boundary conditions) have been analyzed in 
some detail. This report continues in the same spirit. 

2 . Status 

The principal results of the last period are summarized 
as follows: 

a) The complete Karman-Rostovstev plate problem 
was stretched into the correct affine space 
and appropriate solutions for uniform static 
load and for ultimate uniaxial load were pre- 
sented to show the sweeping generality of the 
results . 

b) New material parameters R and X were intro- 

duced, which are combinations of the four 
orthotropic constants; E^, E^ , ^ 12 ' 

These new parameters embody the characteris- 
tics of the fundamental plate parameters D* , 
e and H* . It was found that many common 
orthotropic materials cluster within the 
rather narrow bands 1 < R < 3 and 5 < X < 15 . 
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3. Progress During Report Period 
Progress in the last reporting period has been made 
principally in the following three areas. 

(a) The following P. D. E. describes the behavior of regu- 
lar symmetric cross-ply plates. 
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Examination reveals that this equation differs from the 
specially orthotropic plate equations only in that the con- 
stant, - w , is interchanged with D*. Hence, all solution 
techniques, similarity rules and similarity variables used 
in the orthotropic case are available to facilitate the solu- 
tion of the above cross-ply composite plate equation. 

Few specific solutions exist (all of them in untrans- 
formed real space) for the cross-ply symmetric plate, and 
they are not helpful in understanding how solutions change 
in response to material constant changes. 

The central strong-point of the present formulation, 
the dependence of specific solutions (i.e., static loading, 
vibration, stability) on only one parameter, oi , remedies 
that situation. 

A particularly revealing solution of the vibrating, 
axially-loaded simply supported plate is available in terms 
of the author's similarity variable, T. The results are 
shown in Figure III-E-1 (for all possible axial loadings 
less than the buckling load) for fimn/n^ versus T with a 
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parametric dependence on the parameter ot. Notice that the 
entire frequency spectrum is displayed. 

(b) A chance reading of a 1982 Winter ASME paper (authored 
by S. Batdorf) prompted the author to write an extended note 
using affine transformation theory to make Batdorf ' s equa- 
tions independent of any constants present in his original 
equations. This procedure removes the only known barrier 

to the effective use of Batdorf 's results in assessing local 
stress distributions in damaged composites. 

(c) The rotational derivatives (about some unit vector n) of 
arbitrary order tensors have been derived and detailed re- 
sults have been worked out for fourth order systems, 

In particular, the extremum paths of the elasticity tensor, 

E. ,, can be calculated, and in fact, these derivative re- 
ijkl 

suits are linearly related to the tensor components them- 
selves. This work was motivated by the "derivative-results" 
due to Tsai (and independently derived by the author a few 
months later) ; it seemed natural to expect a generalization 
of these very useful results. 


4. Plans for Upcoming Period 
The author, having taken an extended leave of absence 
from RPI, regretfully must end his association with the 
NASA/AFOSR research effort. This association afforded the 
author ideal conditions for creative work, and the ever- 
steady support of the NASA and AFOSR project monitors is 
gratefully acknowledged. 
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"The Affine Equivalence of Local Stress Distributions in 
Damaged Composites and Batdorf's Electric Analogue" 

To be published in the AIAA Journal, July 1983. 

"Generic Buckling Curves for Specially Orthotropic Rectanular 
Plates", with G. Oyibo. 

To be published in the AIAA Journal, August, 1983. 
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PART IV 

GENERIC STRUCTURAL ELEMENTS 


IV- A COMPACT LUG DESIGN 
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IV-A COMPACT LUG DESIGN 

Senior Investigator: D. B. Goetschel 

1. Introduction 

A critical aspect of many aerospace structural elements 
is the load transfer that takes place between the connecting 
lugs at the ends of the structure and the portions of the 
structure wherein the loads are well-distributed if not uni- 
form. Such lugs are often highly loaded and have very com- 
plex stress states. Further, since they must mate with con- 
necting parts, these lugs are usually designed within rather 
stringent dimensional envelope constraints. As compared to 
lugs which are made from high-strength steel , for example , 
meeting the geometric constraints , even with a composite 
structural design making maximum use of unidirectional gra- 
phite-epoxy, has proven to be a difficult task. This re- 
search is intended to solve the load-volume problem in favor 
of composites rather than either reverting to designs using 
200 ksi steel lugs attached to a composite strut or forcing 
the dimensional constraints to be relaxed with a redesign of 
the mating parts . 

The drag strut of the Lockheed L-1011 is a specific 
primary structure which has been taken as an example to in- 
vestigate heavily loaded pinned connections. While it is 
readily assumed that the column aspect of the structure (see 
Figure IV-A-la) could be redesigned of graphite/epoxy, it is 



Figure IV-A-la. Lockheed L-1011 Engine Drag Strut 




DIMENSIONS IN INCHES 


MAXIMUM COMPRESSIVE LOAD 
MAXIMUM LATERAL LOAD - 


(Schematic) 


207.7 KIPS 
- 236.8 KIPS 
12.0 KIPS 


Figure IV-A-lb. Lug Experiment: Full-Scale L-1011 Drag Strut 
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not. so clear as to whether a graphite/epoxy pinned lug con- 
nection can be designed to withstand the design loads, given 
the original geometric constraints. 

2 . Status 

A number of lug configurations using ply layups (see 
Figure IV- A- lb, for example) and several means intended to 
force failure modes into net tension were fabricated and 
tested. Specimens included quasi-isotropic strips [0/±45/90] 

s 

not scale:d to any particular design, [0/±45] quarter-scaled 

s 

to the L-1011 engine drag strut lug; with round and rectangu- 
lar ends; with and without stress relieving, longitudinal 
slots at the ends of the holes; with and without steel bush- 
ings in the bolt-loaded holes; with and without cap strips. 
Twelve tests specifically comparing slotted and unslotted 
lugs all resulted in failure in combined shear-out and net- 
tension. Although the failure loads all fell within roughly 
the same scatter, all specimens without slots showed a more 
nearly pure tension failure. This showed that the slots im- 
prove stiress distributions in the direction intended. 

3. Progress During Report Period 
a. Stress Relief Cut . 

The stress relief cut concept which was applied to the 
quarter-scale lug (72-ply) models and reported in the previous 

ni* 

progress report was further investigated. Since, in the 

* 

Numbers in brackets in this section refer to the references 
which are listed on page 109. 
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earlier work both the slotted and unslotted (see Figure IV- 
A-2a and 2b) lugs exhibited failure inodes involving shear-out 
at similar load levels , the stress redistribution effect ex- 
pected from the slotted geometry was not verified. The oc- 
currence of a failure more closely approaching a net tensile 
mode would show the benefit of the slots, if such benefits 
exist . 

To isolate the effect of the slots, twelve 16-ply 
(quasi-isotropic layup [0/45/-45/90] 2 g) specimens with rec- 
tangular extended heads were tested during the current per- 
iod (see Figure IV-A-2c) . The squared ends were intended to 
force tensile failure. Bearing failure unexpectedly occurred 
at 955 pounds, however, with the first four specimens tested, 
a 35% decrease in bearing strength is shown to be expected 
for these 16-ply specimens as compared to 72-ply specimens, 

biT 

according to F vs. D/t plots given in the advanced compos- 

(21 

ites design guide . Thus , the premature bearing failures 
can probably be attributed to increasing diameter to thick- 
ness ratio. 

[3] 

Similar testing done, at the University of Delaware 
has shown that bearing stresses, in similar specimens, are 
much higher than those predicted by finite element modeling. 
Three-dimensional effects are suspected as playing a large 
role in this phenomenon. This appears even more likely in 
the case of thin laminates, where free edge stress may cause 
break-away of the outer plies, reducing the effective bearing 
surface by a greater percentage than for thick specimens. 
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To overcome the D/t effect, the area of the net tension 
plane was reduced by reducing the width from 1 . 3 inches to 
0.7 inches on the remaining eight specimens. Net tension 
failure then was clearly evident in each of the subsequent 
tests. The slotted lug specimens showed, on the average, 

21% higher ultimate load capability than unslotted specimens. 
Further, there is no overlap in the data scatter (see Table 
IV-A-1) . This result strongly suggests that the stress re- 
lief cuts are effective in reducing tensile stress concen- 
trations at the edge of the loaded hole. It should be noted 
that the final width to diameter ratio (W/D = 1.4) is shown 

to result in relatively lower tensile stress concentra- 

[4 51 

tions ' than the original width [Note: Here tensile 

stress concentration is defined as peak tensile stress di- 
vided by av... net tensile stress, = (a0) The 

effect of slots may be even more pronounced in geometries 
with larger W/D ratios. Despite the success of these tests, 
the stress relief cut modification appears to be of rather 
limited application. In terms of increased ultimate strength, 
it appears most effective for pinned connections with large 
W/D that would otherwise still fail in tension. It has been 
shown to have negligible effect on shear failures and is 
vulnerable to transverse loading through the pin. 

Attempts to verify these results with finite element 
analyses were frustrated by the assumption of a cosine normal 
load distribution in the available program, and such is not 
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TABLE IV-A-1 

SUMMARY OF FAILURE INFORMATION; 
SLOTTED AND UNSLOTTED TENSION LUGS 


Specimen Description 


Unslotted 


II II 


3 

4 


II 


II 


It 


II 


5 Slotted 

6 " " 


8 


II 


Failure Load Lbs, 
(psi) 


359.6 (40,864) 
311.3 (35,375) 
338 (38,410) 

364.1 (41,375) 

410.2 (46,610) 

429.3 (48,780) 
383.2 (43,550) 
433.8 (49,295) 


Mode of 
Failure 


Net tension 


Av. 

Unslotted 

343 

(38,980) 

Av. 

Slotted 

414 

(47,060) 


All specimens [0/45/-45/90] 2 ^ 



100 


valid for W/D < 2.0^^^. The measured failure stresses were 
roughly 30% lower than the peak stresses predicted by the 
finite element model, probably because of the inaccuracy of 
the model for this particular geometry. 

b. Capstrips 

Another approach pursued in this research was the use of 
a partial capstrip made from unidirectional strips of gra- 
phite/epoxy to strengthen lug designs. The concept is that 
the stiff capstrip will cause a stress redistribution, re- 
ducing the peak stress at the hole edge and raising the aver- 
age stress, at which the inner laminate will fail. 

Our first attempts were reported in the previous pro- 
gress report. In that case, a .1-inch capstrip (see Figure 

IV-A-3a) was bonded to an inner laminate of [0/+45] ^ ^ . 

12s 

The resulting mode of failure involved the capstrip shearing 
off along one side of the lug. Whether or not this precipi- 
tated a subsequent failure, or whether it occurred after 
failure of the inner laminate, was not clear. This specimen 
failed at a load of 12,750 pounds (40,220 psi) . 

It was decided to attempt to wrap the capstrip entirely 
around the inner laminate, forming a tension loop (see Figure 
IV-A-3b) to reduce the tendency of the capstrip to shear off. 
There are, of course, manufacturing difficulties. Whenever 
a laminate is laid up on a surface with curvature and then 
cured, the laminate thickness changes by a significant per- 
centage, as the resin bleeds out, resulting in slack material. 
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Some means must be applied to prevent resin pockets from 
occurring on the curved inner surface , since such imperfec- 
tions can.be highly detrimental. 

A fixture was developed (see Figure IV-A-4) which ef- 
fectively pulls out any slack to the bottom of the lug 
(see Figure IV-A-5a) . This fixture was found to work very 
well in eliminating the slack. It was difficult to tightly 
bag the entire apparatus for cure, however, and as a result, 
uneven pressure on the capstrip caused uneven capstrip 
thickness through the depth of the lug (see Figure IV-A-5b) . 

Despite such fabrication imperfections, tests were con- 
ducted to determine the effect of capstrip width on strength 
and failure mode. These tests involved specimens with .1", 
.2" and .3" depth capstrips. The results (see Table IV-A-2) 
are encouraging, although more work will be required to es- 
tablish trends. All three failures were of the shear-out 
type, with the failure occurring on a shear plane extending 
radially about 25 degrees above the net tension plane. The 
failures appeared to initiate in the inner laminate at the 
hole edge. 

Two concepts were considered in an effort to address 
the apparent initiation of failure at the hole edge. The 
first was the use of a steel bushing that could plastically 
deform to redistribute and reduce the peak stresses. The 
second was the use of some material near the hole which has 
a greater fracture toughness and lower stiffness than gra- 
phite-epoxy to reduce the peak stress. 
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Aluminum Tension Strip 

Peel Ply, Porous Release, Bleeder Cloth 

Graphite Capstrip 

Steel Pin - Fixed to Cure Plate 

Pre- cured Inner Laminate 


Wood Spacer 

Steel Bar - Fixed to Cure Plate 
Compression Springs 


Figure IV-A-4. 


Capstrip Cure Apparatus 



4 



(5a) 

Typical Epoxy Gaps and 
Fiber Slack Resulting From 
Curing Around a Corner 


Figure IV-A-5 . Effects of Capstrip Manufacturing Fixture 
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TABLE IV- A- 2 

SUMMARY OF FAILURE INFORMATION; 
LUGS WITH VARIOUS CAPSTRIPS 


Description 


.1" Straight* 
graphite capstrip 

.1" Continuous 
graphite capstrip 

.2" Continuous 
graphite capstrip 

.3" Continuous 
graphite capstrip 

Hybrid capstrip 
with bushing 


Failure Load Lbs. 
(psi) 

12,750 (40,220) 

10,783 (38,238) 

11,512 (40,823) 

10,536 (37,361) 

11,620 (41,205) 


Mode of Failure 


Capstrip shear off/ 


Shear 

out 

at 

hole 

Shear 

out 

at 

hole 

Shear 

out 

at 

hole 

Shear 

out 

at 

hole 

Splitting 

of 

cap- 


strip 


* 

This specimen [0/±45]^2g' 

All other specimens [0„/±45]_ 

Z oS 
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A combination of these concepts was tried in a hybrid 
capstrip specimen (see Figure IV-A-6) , incorporating a 
1/16” steel bushing with a .169" inner capstrip of unidi- 
rectional fiber glass/epoxy and a .169" outer capstrip of 
unidirectional graphite-epoxy. To minimize the capstrip 
slack around the curved end of the lug, the specimen was 
cured in two cycles, one for each successive capstrip. 

An initial failure was reached at 11,620 pounds (41,200 
psi) , and load application halted at that point. Cracks had 
propagated through the capstrip in the plane of the lug. 

The addition of ±45° plies or metallic layers interspersed 
throughout the capstrip would deter this mode of failure and 
result in a higher strength. However, because the individual 
contributions of the bushing, fiber glass and graphite cap- 
strips are undeterminable from one specimen involving all 
three , it can only be said that an optimum lug design would 
probably incorporate some degree of each. 

c. Computational Aspects 

An important attribute of the design procedure is its 
ability to predict failure loads in preliminary stages so as 
to prevent the squandering of manufacturing and testing re- 
sources in the pursuit of invalid designs. The method of 
analysis used extensively in this study is the finite ele- 
ment package, "POFES" (People Oriented Finite Element System) , 
developed at RPI by Dr. M. S. Shephard. The model being used 
in the process of the research has advanced in sophistication 
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Figure IV-A-6 . Hybrid Capstrip with Steel Bushing 
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by stages. Initially, a constant strain triangle mesh with 
a discrete cosine-normal load distribution and isotropic 
material properties was used. Now, second order quadrila- 
teral elements with a continuous, work equivalent, cosine- 
normal load distribution generated by a program developed in 
this study and anisotropic material properties are available 

in the model. Verification of the model with the results of 
f71 

De Jong has been quite successful. 

It is still recognized, however, that neither the re- 
sults of our finite element analysis, nor De Jong's for that 
matter, provide exact replications of the actual stress field 
around the pin-loaded hole. The necessity of assuming a load 
distribution at the bolt-hole interface and of ignoring three- 
dimensional effects certainly can be expected to lead to in- 
exactness. However, it has been shown that use of an assumed 
load distribution such as the cosine normal distribution can 
produce stress results which correlate reasonably with ex- 
perimental failure data when averaged over a certain charac- 

f 8 1 

teristic distance . We are rigorously testing this stress 
averaging technique for application to bolted joints. 

4 . Plans for Upcoming Period 
It appears from our studies to date, that no feasible 
composite lug design is likely to reach high density design 
criteria such as imposed by Lockheed for the L-1011 drag 
strut pinned connection. A graphite lug constrained to the 
same maximum dimensions as a high strength steel lug doesn't 
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seem capable of withstanding as high a tensile load. Com- 
posite pinned connections as typified by the lug studied 
here , however , are potentially advantageous under more gen- 
eral design envelopes. Despite the lower load/volume effi- 
ciency of graphite/epoxy in this case, it is clearly more 
efficient in terms of load/weight. Our research indicates 
that the graphite lug is nominally 85% as strong as the 200 
ksi steel lug, yet it only weighs 20% as much as the same 
steel lug. The composite thus demonstrates a load/weight 
efficiency 425% greater than steel. Thej weight savings re- 
sulting from use of graphite lugs would seem to justify their 
use in aircraft design. 

This project will be concluded with the work reported 
here and further predictive analytical studies making use of 
averaged stress and characteristic distance techniques. 
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PART V 

PROCESSING SCIENCE AND TECHNOLOGY 


V-A INITIAL SAILPLANE PROJECT: 
V-B SECOND SAILPLANE PROJECT: 


RP-1 

RP-2 
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V- A INITIAL SAILPLANE PROJECT; THE RP--1 

Senior Investigators: F. P. Bundy 

R. J. Diefendorf 
H. Hagerup 
H. Scarton 

1 . Status 

Recognizing that the resins in composite materials can 
change properties somewhat with aging and exposure to tem- 
perature/moisture cycles, the wing/fuselage ensemble struc- 
ture of the RP-1 glider has been static tested about once a 
year since its initial fabrication in '79- '80. This initial 
series of tests consisted of simple vertical bending, and 
the last of this series was done on August 5, 1981. That 
test consisted of simple bending using a distributed load 
up to a maximum of 590 pounds per wing, which corresponds 
to a little over 4 g's for average pilot weight. 

Some changes were made in the wing beam "carry- through" 
structure of the RP-2 glider, in the suinmer of 1982, with a 
eye to making it stronger and more fail-proof. In subsequent 
tests of the RP-1 glider, the wing-fuselage assembly was 
subjected not only to simple vertical bending, as in pre- 
vious years but, in addition, to a bending-torsion test with 
fore-aft bending components, which would simulate the 
stresses in a high speed pull-out. It was in this mode that 
the RP-2 wing beam carry-through failed (see Section V-B of 
the 43rd Semiannual Progress Report, dated December 1982) . 
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In the simple bending tests, the fuselage was held in 
a horizontal position and the weights were laid directly 
over the wing beam. In the bending/torsion tests, the fuse- 
lage was positioned on a ten degree slant (nose low, tail 
high, inverted) and the weights placed on a line fourteen 
centimeters aft of the wing beam. Thus, both torsion and 
fore-aft bending were induced, in addition to vertical bend- 
ing. These two static tests were survived very successfully 

Following the static tests, in July of 1982, winch- 
launched flight tests were carried out to evaluate the in- 
fluence of aerodynamic fairings installed earlier that year, 
on sink rate and glide slope. The effect was to lower the 
minimum sink rate to about 200 FPM (from 250 FPM in the un- 
faired configuration) and to increase the optimum glide 
slope from about 10 or 11:1 to 14:1. 

Tests have been conducted on structural elements of two 
different kinds, typical of those used on the RP-1 and 
hidden from view, in the falls of '81 and '82. Five of each 
kind were tested each year. It appears that degradation 
with time is well within the spread in strength due to fa- 
brication or materials quality. 

2. Progress During Report Period and 
Plans for Upcoming Period 

With the achievement during the summer of 1982 of sus- 
tained thermal soaring of the RP-1 on two separate occasions 
the first all-composite glider built under this program has 



115 


been set aside. Some additional flights may be undertaken 
in the future for demonstration purposes , to stimulate stu- 
dent interest in the program, to maintain an experienced 
flight crew, to test ideas for further aerodynamic refine- 
ment and for their influence on structural life. In addi- 
tion, the program of yearly structural proof-tests and 
structural subassembly testa will be continued, for the pur- 
pose of securing information on possible long-term degrada- 
tion of the structure. The emphasis of the work, however, 
has shifted to the second aircraft, the RP-2. 



lie 
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V-B SECOND SAILPLANE PROJECT; THE RP-2 

Senior Investigatiors : F, P. Bundy 

R. J. Diefendorf 
H. Hagerup 
H. Scarton 

1. Status 

The bending- torsion failure and secondary damage exper- 
ienced during the static test in the fall of 1981 necessi- 
tated structural redesign of the center wing-connection and 
extensive^ development of advanced composite material repair 
techniques on our part. Subsequently, the center wing- 
connection was partly rebuilt, and the two damaged outboard 
portions of the wings were completely rebuilt. This involved 
developing a splicing technique in which autoclave-cured 
primary structures (i.e., spar caps and splice material) were 
joined by cold bonding. The fuselage framework and carry-through 
structure; were also modified, both to provide an improved 
load-path and to accommodate the somewhat larger redesigned 
wing-connection. The former modification included spreading 
the pins which transfer wing root rib loads to the fuselage; 
the greater distance between, them reduces the magnitude of 
torsion cmd fore-aft bending induced forces. 

2. Progress During Report Period 
The rebuilt outboard portions of the wings were joined 


with the undamaged inboard portions, and wing skins to cover 
the area, of these repairs were installed and skin splices 
applied. Aileron and flap connections, linkages and control 
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rods were installed. Modifications were made to the root 
ribs, as necessitated by the new shear-pin locations and 
changes in the carry-through structure generally. Rebuild- 
ing of the center wing-connection was completed, including 
fabrication and installation of the heavy U-channels which 
are now an important part of the design. Final adjustments 
were made to the fuselage framework, the wings were con- 
nected to each other and to the fuselage, and the repaired 
airframe was ready for structural proof-tests at the end of 
the reporting period. 

Meanwhile, a wind tunnel model of the RP-2 fuselage, 
including stub wings, was prepared, for the purpose of veri- 
fying that there are no special problems of aerodynamic in- 
teraction at the wing-fuselage intersection. The model is 
ready for wind tunnel testing at this writing. 

In addition to their becoming familiar with advanced 
composite design and fabrication techniques, and working on 
the RP-2 airframe, participants in the Sailplane Project 
also do some original work on the following; 

a) fabrication technique and structural properties 
of tube connections; 

b) fabrication technique and structural properties 
of double wall tubes with honeycomb between the 
walls; 

c) fabrication of I-beam with sinusoidal web and 

d) design and fabrication of a new concept for a 
variable geometry wing-airfoil section. 



119 


3 . Plans for Upcoming Period 
The completed RP-2 airframe will be subjected to struc- 
tural proof-tests later this spring. Non-structural detail 
work on the wings will be completed during the summer and 
fall. The wind tunnel test on wing-fuselage interaction 
will also be performed in late spring, provided that the 
wind tunnel test results are positive, fabrication of the 
female molds for the fuselage will start immediately. The 
fuselage shell will be laid up in early fall, and the final 
detail work on the entire aircraft will have begun by the 
end of the next reporting period. 

4. Current Publications or Presentations by 
Professors Bundy and Scar ton on this Subject 

"Rensselaer's 'home-built' RP-2", F. P. Bundy v/ith S. J. 
Winckler. 

To be published in Soaring, June 1983. 

"Acoustic Emission Proof Testing of Composite Sailplanes", 

H. A. Scarton. 

To be presented at the First International Symposium 
on Acoustic Emission from Reinforced Plastics, CARP, 

San Francisco, CA, July 19-21, 1983. 
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PART VI 

TECHNICAL INTERCHANGE 
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INTRODUCTION 


Technical meetings, on- and off-campus, provide impor- 
tant opportunities for interchange of technical information. 
Because of the large number of composites meetings, a cen- 
tral catalog with all upcoming meetings is being maintained 
and distributed periodically. In this way we help to assure 
that a Rensselaer faculty/staff member can participate in 
important meetings. The calendar for this reporting period 
is shown in Table VI-1. Meetings attended by RPI composites 
program f aculty/staf f/students during the reporting period 
are shown in Table VI-2. Some meetings particularly rele- 
vant to composites, held on-campus with special speakers, 
are listed in Table VI-3. A list of composite-related vis- 
its to relevant organizations by RPI f aculty/staf f/students , 
with the purpose of each visit outlined, is presented in 
Table VI- 4. 

For the last several years, as the diversity of the 
research conducted within this program has increased, once- 
a-week luncheon programs have been held to insure informa- 
tion transfer among the faculty and graduate students in- 
volved (listed in Part VII - Personnel of this report) . 

These meetings are held continuously when classes are in 
session and are known as "Brown Bag Lunches" (BBL' s) , since 
attendees bring their own. Each BBL is an opportunity for 
graduate students and faculty to present briefly plans for. 



TABLE VI-1 


1982 
10/4-7 

10/5-7 

10/7-8 

10/17-21 

10/25-27 

10/25-28 

11/1-4 

11/17 

12/15-17 

12/15-17 

12 / 20-22 

1983 
1/10-12 

1/11-13 


CALENDAR OF COMPOSITES -RELATED MEETINGS 
(September 30, 1982 through April 30, 1983) 


Symposium on Advanced Trends in Structural and 
Solid Mechanics, Washington, D.C. Sponsored by 
George Washington University and NASA Langley. 

8th Annual Mechanics of Composites Review, Dayton, 
OH. Sponsored by W.P.A.F.B. 

High Performance Composites, CFC 8, Toronto, 
Canada. Sponsored by the Canadian Commission for 
Research on the Strength and Fracture of Materials 

14th Annual Tri-Service Manufacturing Technology 
Conference, Phoenix, AZ. Sponsored by the Air 
Force and DOD. 

The Plastics ABC's: Polymer Alloys, Blends and 

Composites, Bal Harbour, FL. Sponsored by SPE 
and NATEC. 

4th International Conference on Composite Material 
(ICCM-IV) , Tokyo, Japan. Sponsored by AIME and 
the Japanese Society for Composite Materials. 

Japanese Carbon Conference, Toyohashi, Japan. 

2nd Chautauqua Conference on Productivity in Engi- 
neering and Design, The CAD Revolution, Kiawal 
Island, SC. 

International Conference on Biaxial/Multiaxial 
Fatigue, San Francisco, CA. Sponsored by ASTM. 

Indian Carbon Conference, New Delhi, India. 

Silver Jubilee Meeting: Recent Advances in Mate- 

rials Research, Bombay, India. 


21st Aerospace Sciences Meeting, Las Vegas, NV. 
Sponsored by AIAA. 

Composites in Manufacturing II Conference and 
Exhibits, Anaheim, CA. Sponsored by SME . 
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1983 

1/16-20 

1/17-20 

1/21 

1/24-27 

2/7-11 

3/18 

3/21-25 

3/23-25 

3/28 

4/10-15 

4/12-14 

4/12-14 

4/20-22 

4/20-22 


•1 continued 


7th Annual Confeirence on Composites and Advanced 
Materials, Cocoa Beach, FL. 

Gordon Conference on Composite Materials, Santa 
Barbara, CA. 

Society for Advancement of Material and Process 
Engineering (SAMPE) Meeting on Composites, San 
Jose, CA. 

6th Conference on Fibrous Composites in Structural 
Design, New Orleans, LA. Sponsored by the Army, 
Navy, Air Force, FAA and NASA. 

38th Annual Conference, Houston, TX. Sponsored 
by the Society of Plastics Industry. 

Applied Mechanics Colloquia, University College, 
Swansea, U.K. 

Meeting of the High Polymer Division of the Ameri- 
can Physical Society, Los Angeles, CA. 

National Composite Structures Specialists Meeting, 
Philadelphia, PA. Sponsored by AHS. 

American Society Of Metals Symposium, Utica, NY. 

56th Panel Meeting/Specialists Meetings, Structures 
and Materials Panel of AGARD, London, U.K. Spon- 
sored by NATO. 

Fiber Producers Symposixam and Exhibition, Green- 
ville, SC. 

28th National Symposixam, Anaheim, CA. Sponsored 
by SAMPE. 

Symposixam on Computer-Aided Geometry ModelING, 
Hampton, VA. Sponsored by NASA Langley. 

National and International Meeting on Composite 
Materials, Milan, Italy. Sponsored by AMME. 
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TABLE VI -2 

COf/LPOSITES-RELATED TECHNICAL MEETINGS ATTENDED OFF-CAMPUS 
(September 30, 1982 through April 30, 1983) 

1982 

10/4-7 Symposium on Advanced Trends in Structural and 
Solid Mechanics (Prof. Loewy and Grad. Student 
F. Luban-Tobing) , Washington, D.C. 

Professor Loewy chaired session on Compos- 
ite Materials and Structures. 

F. Liimban-Tobing presented the paper: 

"Finite Element Modeling of Moisture Ef- 
fects in Graphite/Epoxy Composites". 

10/4-7 8th Annual Mechanics of Composites Review (Prof. 
Goetschel) , Dayton, OH. 

10/25-28 4th International Conference on Composite Mate- 
rials (ICCM-IV) (Prof. Diefendorf) , Tokyo, Japan. 

Professor Diefendorf chaired a session and 
presented the paper, "Residual Stress in 
High Modulus Carbon Fibers". 

Professor Hoff presented the paper , 

"Stress Concentrations in Cylindrically 
Orthotropic Plates with Radial Varia- 
tion of the Compliances" (Co-authored 
with C. Muser) . 

11/1-4 Japanese Carbon Conference (Prof. Diefendorf), 
Toyohashi, Japan. 

Professor Diefendorf presented the lec- 
tures : 

"Mesophase Formation in Polynuclear 
Aromatic Compounds" 

"Mesophase Formation from Fractions 
of a Mesophase Pitch" 

"Molecular Weight and Molecular Weight 
Distributions in Pitch" 

"The Chemical Vapor Deposition of Car- 
bon on Graphite Surfaces" 


"k 

Work supported by the subject grant. 
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TA BLE VI" 

1982 

11/17 

12/15-17 

12/15-17 


12 / 20-22 

19 83 
1/17-20 

1/21 

3/18 


2 continued 


2nd Chautauqua Conference on Productivity in Engi- 
neering and Design, The CAD Revolution (Prof. 
Shephard), Kiawal, SC. 

Professor Shephard presented the paper, 

"The Effective Use of Finite Element 
Methods in Engineering Design" . 

International Conference on Biaxial/Multiaxial 
Fatigue (Prof. Krempl) , San Francisco, CA. 

Professor Krempl presented the paper, 

"Biaxial In-Phase and Out-of -Phase Be- 
havior of Graphite-Epoxy Tubes". 

Indian Carbon Conference (Prof. Diefendorf) , New 
Delhi, India. 

Professor Diefendorf presented two papers 
and a plenary lecture as follows : 

"Transverse Thermal Expansion of Carbon 
Fiber/Epoxy Matrix Composites" 

"Residual Stresses in Carbon Fibers" 

"Mesophase Pitches" 

Silver Jubilee Meeting; Recent Advances in Mate- 
rials Research (Prof. Diefendorf), Bombay, India. 

Professor Diefendorf chaired a session 
and presented the paper, "Carbon Fibers, 

Present Properties and Future Applica- 
tions" . 


Gordon Conference on Composite Materials (Profs. 
Diefendorf and Sternstein) , Santa Barbara, CA. 

Society for Advancement of Material and Process 
Engineering (SAI4PE) Meeting on Composites (Prof. 
Diefendorf), San Jose, CA. 

Professor Diefendorf presented the Keynote 
Lecture, "RP-X Composite Sailplanes". 

Applied Mechanics Colloquia, University College 
(Prof. Shephard) , Swansea, U.K. 

Professor Shephard presented the paper, 
"Developments in Mesh Generation and Adap- 
tive Techniques". 
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TABLE VI- 

198 3 

3/21-25 

3/23-25 

3/28 

4/10-15 

4/12-14 


2 continued 


Meeting of the High Polymer Division of the Ameri- 
can Physical Society (Prof. Sternstein) , Los Angeles, 
CA. 

Professor Sternstein presented the paper, 
"Characterization of High Performance 
Composites". 

AHS National Composite Structures Specialists 
Meeting (Profs. Diefendorf, Goetschel and Loewy, 
and Mr. Paedelt) , Philadelphia, PA. 

Professor Loewy chaired session, "VI: Anal- 

ysis II". 

American Society of Metals Symposium (Prof. Diefen- 
dorf) , Utica, NY. 

Professor Diefendorf presented, "Princi- 
ples of Composite Materials". 

56th Panel Meeting/Specialists Meetings, Structures 
and Materials Panel of AGARD (Prof. Loewy) , London, 

U.K. 

Fiber Producers Symposium and Exhibition (Prof. 
Diefendorf) , Greenville, SC. 

Professor Diefendorf presented the 
Invited Speech, "A Comparison of High 
Modulus Fibers " . 
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TABLE VI -3 

COMPOSITES-RELATED MEETINGS/TALKS HELD AT RPI 
(September 30, 1982 through April 30, 1983) 


Topic Date Speaker (s) 

1982 

Computer Aided Mechanism 10/12 Arthur G. Erdman 

Professor, University of 
Minnesota 

G. I. Schueller 
Institut fiir Mechanik 
Universitat Innsbruck 
Austria 


Some Considerations of 10/18 

Structural Reliability 
with Respect to Risk Anal- 
ysis of Large Technological 
Systems 


1983 

Design, Manufacturing and 1/26 
Testing of High Speed Ro- 
tating Graphite/Epoxy 
Shafts , 

Biaxial Fatigue of Graph- 2/3 
ite/Epoxy [±45]^ Tubes 


Structural Analysis of 3/1 

Jet Engine Combuster 

Liners 

Recent Trends in the Com- 4/8 
puter-Aided Control of 
Material Processing 
Methods 

Commemorative Sympositim: 4/5 

Rensselaer Salutes 200 
Years of Manned Flight 

daVinci's VTOL Comes 
of Age 


A/C Development in the 
Skunk Works 


Olivier A. Bauchau 
St. Gobain Recherche 
Paris, France 


Tyan-Min Niu 

Doctoral thesis dissertation 
Department of Mechanical 
Engineering, Aeronautical 
Engineering and Mechanics 

Kevin P. Walker 
Habbitt Karlsson and 
Sorensen, Inc. 

Daeyong Lee 

General Electric Corporate 
Research and Development 
Schenectady, New York 


Robert Zincone 
Senior, Vice President 
for Engineering 
Sikorsky Div., UTC 

Ben Rich 

Vice President and General 
Manager 

Lockheed California 
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TABLE VI- 3 continued 


Topic 


Twenty-Five Years on 
"This New Ocean" 


The Year of the Shuttle: 
Fifth Annual Hemke Lec- 
ture 

The Effect of Stress Re- 
laxation on the Tearing of 
Elastomers 

Fibrous Skins 


Date Speaker (s) 

1983 

John F. Yardley 
President 

McDonnell Douglas Astro- 
nautics 

Col. Robert Overmeyer 
STS-5 Pilot 
NASA Astronaut 

4/12 R. S. Rivlin 

Lehigh University 


4/26 S. Wainwright 

Duke University 
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TABLE VI~4 

COMPOSITES- RELATED VISITS TO RELEVANT ORGANIZATIONS 
by RPI Faculty /Staff/Students 
(September 30, 1982 through April 30, 1983) 


Visited 

Date 

By 

Purpose 


1982 



Techniche Hocho- 
chules , Darmstat 
and Braunschweig 
Federal Republic 
of Germany 

10/18, 

20 and 

21 

M. V. Paedelt 

Discussed operations 
of "Akaflieg" groups 
on these campuses 

Toray , Ind . 
Shiga, Japan 

10/26 

Prof. R. J. 
Diefendorf 

Plant visit to one of 
the largest manufac- 
turers of carbon fi- 
bers 

NASA Lewis Re- 
search Center 
Cleveland, OH 

11/19 

Prof. R. J. 
Diefendorf 

Participated in SiC 
Semiconductor Work- 
shop 

Air Force Insti- 
tute of Technol- 
ogy , Dayton , OH 

12/28 

Prof. E. J, 
Brunelle 

Gave invited lecture: 
New Research Areas in 
Composite Structures 


1983 



Sikorsky Divi- 
sion, United 
Technologies 
Corp, Stratford, 
CT 

1/27 

Profs . R. J . 
Diefendorf 
and D. Goet- 
schel 

Made plant tours and 
discussed potential 
composite research 
projects 

NASA Langley Re- 
search Center 
Hampton, VA 

1/31- 

2/4 

Prof. M. S. 
Shephard 

Attended workshop on 
finite element model- 
ing and presented 
talk: "Finite Element 
Mesh Generating Tech- . 
niques and the Inter- 
face of Finite Elements 
and Geometric Modeling" 

University of 

Michigan 

Ann Arbor , MI 

3/14 

Prof. S. S. 
Sternstein 

Presented talk : 
"Characterization of 
Composites and Neat 
Resins" 



problems encountered in and recent results from their indi- 
vidual projects. These sessions also are occasions for 
brief administrative reports , usually on the part of one of 
the Co-Principal Investigators, and for brief reports on 
what transpired at off-campus meetings attended by one of 
the participants (as listed in Tables VI-2 and VI-4 of this 
report). Visitors from off, campus, who are at RPI during 
a BBL day are often invited to "sit in". A calendar of 
internal, oral progress reports as they were given during 
the last reporting period at BBL's is listed in Table VI-5. 
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TABLE VI- 5 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
BROWN. BAG LUNCH (BBL) SCHEDULE 
(September 30, 1982 through April 30, 1983) 


Responsible 


Date 

Topic 


Faculty 

1982 




10/1 

Administrative Report 

R. 

Loewy 


Computer-Aided Design and Analysis 

M. 

Shephard 


Resin Matrix Research 

J. 

Dief endorf 

10/8 

Administrative Report 

J. 

Dief endorf 


Generic Structural Components 

D. 

Goetschel 


Fcitigue Studies Research 

E. 

Krempl 

10/15 

Administrative Report 

J. 

Dief endorf 


Fcibrication Technology Experiments 

F. 

Bundy 



H. 

Hagerup 



V. 

Paedelt 


Report on Mechanics of Composites Review 

D. 

Goetschel 

10/22 

Aclministrative Report 

R. 

LoevTy 


Matrix Characterization and Environ- 
tcil Effects 

S. 

Sternstein 


Computer-Aided Design and Analysis 

M. 

Shephard 

10/29 

Administrative Report 

R. 

Loewy 


Generic Structural Components 

D. 

Goetschel 


Advanced Analytical Methods 

E, 

Brunelle 

11/5 

Aclministrative Report 

J. 

Diefendorf 


Fabrication Technology Experiments 

F. 

Bundy 



H. 

Hagerup 



V. 

Paedelt 


Fatigue Studies Research 

E. 

Krempl 

1983 




1/28 

Administrative Report 

R. 

Loewy 


Report on Gordon Conference 

J. 

Diefendorf 
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TABLE VI-5 continued 


Responsible 


Date 

Topic 

Faculty 

1983 



2/4 

Matrix Characterization of Environ- 
mental Effects 

S. Sternstein 


Fabrication Technology Experiments 

F . Bundy 
H . Hagerup 
V. Paedelt 

2/11 

Generic Structural Components 

D. Goetschel 


Computer-Aided Design and Analysis 

M. Shephard 

2/18 

Resin-Fiber Interface Research 

J. Diefendorf 


Composites Fatigue Research 

E . Krempl 

2/25 

Administrative Report 

R. Loewy 


Pin-Loaded Holes in Composite Plates 

W. Kim 


Computer-Aided Design and Analysis 

M. Shephard 

3/4 

Matrix Characterization and Environ- 
mental Effects 

S. Sternstein 


Fabrication Technology Experiments 

F . Bundy 
H. Hagerup 
V. Paedelt 

3/11 

Generic Structural Components 

D. Goetschel 


Advanced Analytical Methods 

E. Brunei le 

3/25 

Report on Vertol Research 

E. Krempl 
W. Kim 


Resin-Fiber Interface Research 

J. Diefendorf 

4/1 

Administrative Report 

S . Wiberley 


Composites Fatigue Research 

E. Krempl 


Computer-Aided Design and Analysis 

M. Shephard 

4/8 

Matrix Characterization and Environ- 
mental Effects 

S. Sternstein 


Fabrication Technology Experiments 

F. Bundy 
H. Hagerup 
V. Paedelt 
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TABLE VI-5 continued 


Date Topic 

19 83 

4/15 Administrative Report 

Generic Structural Components 
Free Edge Failures of Laminates 

4/22 Report on AGARD Meeting 

Resin-Fiber Interface Research 

4/29 Administrative Report 

Composites Fatigue Research 
Computer-Aided Design and Analysis 


Responsible 

Faculty 

J. Diefendorf 

D. Goetschel 
S . Sham 

R . Loewy 
J. Diefendorf 

R . Loewy 

E . Krempl 
M. Shephard 
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PERSONNEL 


Co-Principal Investigators 

Ansell, George S., Ph.D. 
Loewy, Robert G. , Ph.D. 
Wiberley, Stephen E. , Ph.D. 

Senior Investigators 

Brunelle, E. J. , Jr., Sc.D. 
(Aeroelastic and structural 
design and analysis, applied 
raechanics of composite 
structures) * 

Bundy, F. P. , Ph.D. 

(Physical chemistry and 
structures testing) * 

Diefendorf ^ , R. J. , Ph.D. 
(Fcibrication , resin matrix, 
fiber behavior, interfaces)* 

Feeser^. L. J. , Ph.D. 
(Computerr applications and 
graphics , computer-aided 
design, optimization)* 

Goetschel, D. B. , Ph.D. 
(Structural analysis design 
and testing) * 

Hagerup, H. J. , Ph.D. 
(Aerodynamics , configura- 
tion, pilot accomodation, 
flight testing) * 

Krempl , E . , Dr . Ing . 

(Fatigue studies , failure 
criteria) * 

Scarton, H. , Ph.D. 

(Acoustic emission NDE) * 


Dean, School of Engineering 
Institute Professor 
Professor of Chemistry 


Associate Professor of 
Aeronautical Engineering 


Research Professor of Mate- 
rials Engineering 

Professor of Materials 
Engineering 

Professor of Civil Engineering 


Assistant Professor of 
Mechanical Engineering 

Associate Professor of 
Aeronautical Engineering 


Professor of Mechanics and 
Director of Cyclic Strain 
Laboratory 

Associate Professor of Me- 
chanical Engineering and Me- 
chanics 


Fields of Speciality 

^Member of Budget Committee together with Co-Principal Inves- 
tigators 
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Senior Investigators 

Shephard, M. S., Ph.D. 
(Computer graphics, finite 
element methods) * 

Sternstein^ , S. S., Ph.D. 
(Failure analysis, matrix 
behavior, moisture effects)* 


Associate Director, Center for 
Interactive Computer Graphics 
and Assistant Professor of 
Civil Engineering 

William Weightman Walker 
Professor of Polymer Engineer- 
ing 


Research Staff 

Manager & Master Technician, Composites Laboratory 
Paedelt, Volker 
Research Associate 
Kim, Wonsub, Ph.D. 


Graduate Assistants 
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